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ABSTRACT 

We present the modeling of SINFONI integral field dynamics of 18 star forming galaxies at z ~ 2 
from Ha line emission. The galaxies are selected from the larger sample of the SINS survey, based on 
the prominence of ordered rotational motions with respect to more complex merger induced dynamics. 
The quality of the data allows us to carefully select systems with kinematics dominated by rotation, 
and to model the gas dynamics across the whole galaxy using suitable exponential disk models. We 
obtain a good correlation between the dynamical mass and the stellar mass, finding that large gas 
fractions [Mgas ~ Mt,) are required to explain the difference between the two quantities. We use the 
derived stellar mass and maximum rotational velocity Vmax from the modeling to construct for the 
first time the stellar mass Tully-Fisher relation at z ~ 2.2. The relation obtained shows a slope similar 
to what is observed at lower redshift, but we detect an evolution of the zero point. We find that at 
z ~ 2.2 there is an offset in log{M^:) for a given rotational velocity of 0.41 ± 0.11 with respect to the 
local Universe. This result is consistent with the predictions of the latest N-body/hydrodynamical 
simulations of disk formation and evolution, which invoke gas accretion onto the forming disk in 
filaments and cooling flows. This scenario is in agreement with other dynamical evidence from SINS, 
where gas accretion from the halo is required to reproduce the observed properties of a large fraction 
of the z ~ 2 galaxies. 
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1. INTRODUCTION 

In recent years our observational knowledge of the high 
redshift Universe has increased enormously. A variety 
of selection techniques have been developed to create 
well-defined samples of z=l-6 galaxies, from UV to sub- 
millimet er wav elengths (e.g. Gia valisc o et al. l2002l Stei- 
del et al. l2004l Adelberger et al. [2004', Franx et al. 1200^ 
Chapman et al. l2005l . Daddi et al. 2004b). These com- 
plementary and partially overlapping samples are now 
covering the critical redshift range around z 2, which 
is the epoch of major star formati on and mass assem- 
bly in galaxies (Hopkins & Beacom 20061 ). Most of the 
first studies of these high-z populations came from pho- 
tometric measurements and comparisons with template 
spectral energy distributions (SEDs). As a consequence, 
little was known about the detailed properties and dy- 
namics of these systems, since such information can best 
be obtained through sensitive, high-resolution near-IR 
(rest frame optical) spectroscopy. In the last few years, 
using increasingly sensitive IR detectors coupled with 
high throughput spectrometers and integral field tech- 
niques, it has become possible to carry out detailed imag- 
ing spectroscopy of the strongest rest-frame optical emis- 
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sion lines, allowing the assembly stage of massive galaxies 
to be studied directly through observations of the dy- 
namical and detailed physical properties of the d ifferent 
popul ations at z = 1.5 — 3 (see e.g. Erb et al. l2006al . 
l2006c). Indeed, the dynamical state of distant galaxies is 
the key property that enables a quantitative assessment 
of the main drivers of mass assembly and evolutionary 
processes. 

In this context we have initiated SINS, the high-z 
galaxy S'pectroscopic Imaging survey in the N ear-IR 
with SINFONI (Forster Schreiber et al. 120081 ). mak- 
ing use of the capabilities of spatially resolved integral 
field spectroscopy to obtain new insights on the inter- 
nal dynamics, sizes and morphology, dynamical masses, 
metal abundances and stellar populations of high-z mas- 
sive star-forming galaxies. With the direct measurement 
of these quantities, the first study of detailed dynamical 
properties of these sources has been successfully carried 
out. We find evidence for large, massive, gas-rich rotat- 
ing disks in place at z ~ 2, requiring ra pid an d efficient 
mass accretion (Forster Schre iber e t al. 120061 : Genzel et 
al. [2006, 2008; Shapiro et al. 120081: see also e.g. Wright 
et al. i2007. . Bournaud et al. l2007l . van Starkenburg et 
al. 120081) . While major (< 3 : 1) mer gers m ay form such 
disks in some cases (Robertson et al. 12006 ; Robertson & 
Bullock 2008) , on average a more continuous and steady 
mass accretion mechanism, such as throug h col d flows 
along cosmi c web filaments (Keres et al. l2005l Dekel 
& B irnbo im [20M . Dekel et al. 200^, Sommer-Larsen et 
al. l2003l ). could be invoked. This leads to high star 
formation rates with longer time scales than the short 
bursts driven by mergers, and possibly induces the for- 
mation of central bulge s via internal/secul ar dyn amical 
processes ( e.g. N oguchi IT999I . Immeli et al. l2004 Bour- 
naud et al. l2007f ). The latest generation of cosmological 
semi-analytical and hydrodynamical simulations in the 
A-CDM paradigm indicate that minor mergers and/or 
smooth accretion are likely to dominate the assembly of 
the high-z populations, while major mergers would play 
a role in the smaller fraction of hyperluminous sources 
(see e.g. Kitz bichle r & White 2007, Naab et al. 2007, 
Guo & White l2008l Genel et al. 120081 ). A comparison 
of the dynamical properties of different sa mples o f high-z 
galaxies observed in SINS (Bouche et al. l2007| ). shows 
that optically and IR-selected galaxies up to z = 2.2 
follow a similar velocity-size relation as z = late-type 
spirals, in strong contrast with sub- mm selected galaxies 
(see Tacconi et al. i2008i) . 

The spatially resolved integral field kinematics of high- 
z galaxies has therefore proven to be a powerful tool to 
study the major processes connected to the formation 
and evolution of galaxies in the ea rly universe (see also 
Wright et al. ,2007, . Law et a l l2007l Bournaud et al. l2008l 
and van Starkenburg et al. 12008^ . A robust measure of 
the internal kinematics of high-z galaxies is therefore a 
fundamental step in our understanding of the formation 
and early evolution of galaxies. 

In this paper we present quantitative kinematic mod- 
eling of a sample of galaxies at z ~ 2, selected from the 
SINS sample to be the most disk-like according to their 
kinematics, showing a regular velocity "spider-pattern" 
with a clear velocity gradient and the velocity dispersion 
peaking in the center. Thanks to the full two-dimensional 
coverage of our SINS datasets, we are able to model the 



internal dynamics of the warm gas as traced by the Ha 
emission, deriving important constraints on the struc- 
ture of these objects and on the evolution of their scal- 
ing pr operti es, such as the TuUy-Fisher relation (TuUy & 
Fisher il^. 

In Section [2] we describe the sample selection, observa- 
tions and data reduction, along with the extraction of the 
kinematic information from the observed datacubes. In 
Section [3] the disk model used and the fitting technique 
are presented, while the results are discussed in Section|3) 
In Section [5] the stellar mass TuUy Fisher relation at 
z ~ 2 is presented, and our conclusions follow in Sec- 
tion [6l Throughout the paper we assume a A-dominated 
cosmology with Hq — 70 kms~^Mpc^^, flm = 0.3 and 
fl\ — 0.7. In this cosmology, 1" corresponds to ~ 8.2 
kpc aX z — 2.2. All magnitudes are given in the Vega 
photometric system. 

2. SAMPLE AND OBSERVATIONS 

In the framework of the SINS survey (see F orster 
Schreiber et a l [2006I . [2008; Genze l et a l. l200l [20081 : 
Bouche et al. l2007t Shapiro et al. |2008[ ). we used the 
SINFONI near-IR integra l field spectrograph on the E SQ 
VLT (Eisenhauer et al. l2003l . Bonnet et al. l2004( l to 
measure the dynamics of the warm gas in a substantial 
sample of over 60 redshift 1.5 — 3.5 massive star-forming 
galaxies, as traced by their Ha emission redshifted into 
the i7-band (z ~ 1.5) and if-band (z ~ 2.2), or by [QUI] 
at z > 2.8. 

2.1. The SINS high-z disks sample 

Forster Schreiber et al. 2008 describe in detail the 
selection criteria of the SINS z 2 sample. Briefly, 
the SINS galaxies were selected to be massive star form- 
ing galaxies using a variety of complementary techniques 
in UV throu gh IR wavebands (BX/B M criter ion, Adel- 
berger et al. |2004 s -BzK , Daddi et al. I2004U K magni- 
tude, Cimatti et al. l200l IRAC flux, Kurk et al. l2006( l. 
These selection criteria sample fairly luminous {LboI ~ 
IO^-^^Lq) and massive galaxies {M lO^-^-^Mg) 
with star formation rates of 10-300 Mq/jy. These rep- 
resent the bulk of the cosmic star formation activit y and 
stellar mass dens ity at these redshifts ( Reddy et al. l2005t 
Rudnick et al. 120061 : Grazian et al. I2007L Caputi et 
al. l2007t ). Here we present the dynamical modeling 
performed on a subsample of 18 galaxies (see Table (TJ 
Fig. [2]), that are selected for showing ordered rotational 
signatures, i.e. a regular velocity pattern with a clear 
velocity gradient and a velocity dispersion that peaks in 
the center, as well as for their high S/N observations. 

In order to distinguish between rotationally dominated 
galaxies and sources with merger induced complex dy- 
namics, we used, where possible, the method described 
in Sha piro e t al. (|2008[ ). Using kinemetry (Krajnovic 
et al. l2006[ ). we were able to quantify asymmetries in 
both the velocity and velocity dispersion maps of the 
warm gas, in order to empirically differentiate between 
major merging and non-merging (or minor merging) sys- 
tems at high redshift. These criteria take full advantage 
of the wealth of information available with integral fleld 
data, which - unlike broad-band morphology - give di- 
rect information on the dynamical state of the galaxies. 
The method was flrst applied to a subsample of galax- 
ies with the best quality data (see Shapiro et al. l2008f) . 
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and it has been now extended to a a larger subset of 
SINS galaxies (Shapiro et al. in preparation). In total 
11 out of 18 galaxies of the sample discussed here were 
classified as rotationally dominated by kinemetry. For 
the remaining 7 sources (Q1623-BX663, Q2346-BX416, 
D3a-7144, D3a-4751, ZC1101592, GK2471, GK1084), we 
were not able to derive a quantitative classification due 
to the lower S/N of the data. We include them in this 
study since they closely resemble at a visual inspection 
the ones classified as disks in their dynamical properties 
(regular velocity pattern with clear velocity gradient, ve- 
locity dispersion peaking in the center), with no evidence 
of more complex merger induced dynamics. 

The disk like systems selected for this study include 
7 BX/BM galaxies selected by their rest- frame UV col- 
ors dow n to R < 25.5, taken from the Ha survey of Erb 
et al. ()2006d) . 2 galaxies selected with IRAC/Spitzer 
myiB(4.5/iTO) < 23 in the framework of the GMASS sur- 
vey (Kurk et al. l2006f ). and 10 rest frame optical selected 
galaxies meeting the s-BzK color criterion introduced by 
Daddi et al. (2004K )^_These are taken from th e K20 
survey (Daddi et al. l2004aD. De ep-3a (Kong et al. 120061) 
and z-Cosmos (Lilly et al l2007l ) surveys. 

2.2. Observations and data reduction 

We observed t he se lected galaxies with SINFONI 
(Eisenhauer et al. |2003( ) at the VLT UT4 telescope dur- 
ing several campaigns between July 2004 and November 
2007 as part of guaranteed time observations. The H 
or K band grating was used to sample the Ha emission, 
depending on the redshift of the galaxy. Most of the ob- 
servations were carried out in seeing limited mode, using 
the 0.250" X 0.125" pixel scale, which provides a total 
field of view of 8" x 8". The average resolution obtained 
for our 18 galaxies, as sampled by PSF reference star 
observations before and after every hour of target in- 
tegration, is ~ 0.5". The spectral resolution provided 
is about 80 km/s Full Width Half Maximum (FWHM) 
in the K band and 100 km/s in the H band. For one 
source, D3a-15504, the presence of a suitable nearby ref- 
erence star allowed us to use the correction provided by 
Natural Guide Star (NGS) Adaptive Optics (AO) obser- 
vations. For two more objects, ZC782941 and Q2346- 
BX482, the L aser G uide Star (LGS) syst em PARSEC 
(Rabien et al. l2004l . Bonaccini et al. l2006l ) was used to 
provide an average resolution of 0.2" (corresponding to 
only 1.25 kpc at z = 2). In the AO observations of these 
three sources the 0.100" x 0.050" pixel scale was used, 
providing a field of view of 3.2" x 3.2". The observations 
and the data reduction were performed as described in 
Forster Schreiber et al. (I20060 . Abuter et al. (2006), and 
Davies et al. (2007). 

The final result of the data reduction procedures is 
a flux-calibrated datacube, containing the full image of 
the galaxy observed at each wavelength. We are thus 
able to extract from the cube both the one-dimensional 
spectrum at each pixel (or integrated over a larger re- 
gion) and monochromatic images of the field of view at 
different wavelengths. 

2.3. Kinematic extraction 

To increase the S/N ratio, we median filtered the re- 
duced cubes spatially using a FWHM=3 pixels, obtain- 
ing an effective resolution of ^ 0.6", or ~ 0.2" for the 



galaxies observed with the 0.05" pixel scale and AO. 
We then extracted from the smoothed cubes the con- 
tinuum maps, line emission, velocity and velocity dis- 
persion fields of Ha. This was done by fitting a func- 
tion to the continuum-subtracted spectral profile at each 
spatial position in the datacube, masking out the wave- 
length range contaminated with the [Nil] line emission. 
The function fitted was a convolution of a Gaussian with 
a spectrally unresolved emission line profile of a suit- 
able sky line, i.e. a high signal and non-blended line in 
the appropriate band. A minimization was performed 
in which the parameters of the Gaussian were adjusted 
until the convolved profile best matched the data. Dur- 
ing the minimization, pixels in the data that consistently 
deviated more than 3a from the average were rejected, 
and were not used in the analysis. The uncertainties on 
the measured velocity and velocity dispersion maps are 
evaluated through Monte Carlo realizations, perturbing 
the input data assuming Gaussian uncertainties. All de- 
tails and issues of the kinematic extraction with the IDL 
routine LINEFIT are d escribed in more details in Forster 
Schreiber et al. (|2008[ ) and Davies et al. (in preparation). 

3. DYNAMICAL MODELING OF THE GALAXIES 

The observed Ha line emission kinematics of the se- 
lected galaxies can now be studied to robustly quantify 
their dynamical properties. This represents a major step 
forward in the modeling of the dynamics of high-z galax- 
ies, as the full two-dimensional mapping from SINFONI 
is more complete and not biased by a priori assumptions 
about the kinematic major axis and inclination of the 
syste m as in lon g-slit s pectroscopy (see also Wright et al. 
[20f)7l Law et al. l20f)l . 

We compare the velocity and velocity dispersion maps, 
derived from the SINFONI data as described in Sec. [2l 
with an exponential disk model (scale height = 0.1"). 
Although a simple exponential disk could be an over- 
simplification of the structure in some of these high-z 
sources, this represents the best compromise between 
providing homogeneous reliable properties and avoiding 
over-interpretation of the data. A more detailed analysis 
is possible for the higher resolution and S/N datasets, 
using more complex mass distributions (see Genzel et al. 
I2OOI . 

The observed intensity distribution of the Ha emission 
line traces the location of the star forming regions across 
the galaxy, and has a clumpy and irregular distribution 
even for kinematically very regular disks. Therefore the 
line emission light distribution is used in the following 
only to evaluate the size of the emitting region. The best 
fitting disk parameters were derived using an optimized 
minimization routine, comparing the disk model with 
the observed velocity and velocity dispersion extracted 
as described in Sec. 12.31 

3.1. The disk model 

We create the disk models using the IDL code DYS- 
MAL. The code uses a set of input parameters to derive 
a datacube with two spatial and one spectral (velocity) 
axis, from which it is possible to extract 2D morphologi- 
cal and kinematical maps. It begins by creating a face-on 
model of an axisymmetric galaxy with 3 spatial dimen- 
sions [ATo, loi -^o]; in which Zq is normal to the galaxy 
plane. The radial mass profile is specified by the user 
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and can be based on exponential, Gaussian, Moffat, or 
power-law functions. The radial luminosity profile can 
be specified independently, since it can refer to different 
tracers that do not necessarily follow the mass distri- 
bution. The scale height of the disk must also be pro- 
vided. The 3D model is then rotated to the required 
inclination and position angles, to create a cube that 
has axes [Xs,Ys, Zs] where is the projection on 

the sky. Although a warped disk can be modelled by 
specifying these 2 parameters as functions of radius, in 
all cases here the disk was modelled in a single plane. 
To make this coordinate transformation more efficient, 
it is first necessary to identify those pixels in the result- 
ing [Xs,Ys, Zs] cube containing flux, and only fluxes for 
these pixels are calculated. The flux in each pixel is in- 
terpolated based on the sub-pixel location in the original 
[Xo, Yq, Zq] cube. A second matching cube is created in 
which each pixel is assigned an appropriate line-of-sight 
velocity. The velocities are interpolated from a rotation 
curve that is itself generated from the combined mass 
profile (including a black hole mass if one is specified) 
and scaled to the total mass at a given radius. The next 
step is, for each pixel [Xs,Ys, Zs], to generate a Gaussian 
profile to model the line emission, whose scaling depends 
on the flux, whose center depends on the linc-of-sight ve- 
locity, and whose width depends on both the instrumen- 
tal broadening specified by the user and a local isotropic 
velocity dispersion depending on the z-scale height. This 
z- velocity dispersion aoi is estimated in one of two ways: 

1. using an approximation appropriate for large thin 
disks, in which it depends on the scale height and 
mass surface density (Binney & Tremaine i2008t, eq. 
4.302c): 

„„, . fWK ,1, 

where is the scale height and v{R) the rotational 
velocity at radius R; 

2. assuming that the vertical motion is restricted to 
the inner part of an extended mass distribution and 
does not feel the entire gravity: 



which is more suitab le for a thick or compact disk 
(Genzel et al. I2Q08D . and it was used in the mod- 
eling discussed below. 

The total z-vclocity dispersion in the model is then given 
by ctq = ■\/ + (702 , where cto2 is an additional compo- 
nent of isotropic velocity dispersion throughout the disk 
that can be added by the user. This term was intro- 
duced in DYSMAL to account for instrumental broaden- 
ing. It can also be used to include any additional pos- 
sible contribution to the observed velocity widths from, 
e.g., rapid mass accretion, strong feedback from star for- 
mation, or unresolved non-circular motions (hereafter 
globally called "random motions" for simplicity). This 
yields a line-of-sight velocity profile for each pixel in 
[XgjYs, Zs]- Then, at each [Ars,^;] location, the veloc- 
ity profiles for all the [Zs] positions are summed to yield 
the full line-of-sight velocity distribution at that point. 
This creates a cube [Xs, Ys, V] with two spatial and one 



velocity axes. The final step is to convolve each [Xs, Ys] 
plane in this cube with the spatial beam and to rebin to 
the desired pixel scale. 

The reason for beginning with a galaxy in 3 spatial di- 
mensions is to ensure that there is always enough spatial 
resolution along the projected minor axis to yield a re- 
liable line-of-sight velocity distribution. This can be an 
issue when the inclination is high, and in such cases a 
direct 2D projection is insufficient. It also enables one 
to apply a mask that modulates the flux distribution, 
for example to simulate a knotty or clumpy structure. 
The code allows one either to generate a mask contain- 
ing random knots; or to input one if, for example, a more 
ordered pattern such as spiral arms is required. In either 
case, the mask is applied to the galaxy while it is face-on, 
in the [Aro,lo] plane. For this high-z study, however, a 
plain exponential disk was used to keep the number of 
free parameters as low as possible. 

The velocity dispersion across the galaxy disk is an 
important parameter to consider. When there are signif- 
icant random motions, as implied by the large observed 
dispersion, they should also support some of the mass. 
This is currently not implemented in DYSMAL, which 
assumes that all the mass is supported by ordered rota- 
tion in the disk plane. Thus, if v/ao < \/3, the total 
system mass may be underestimated by a factor of 2 or 
more. The average V/ctq for our sample is ^ 4.4, pro- 
ducing an effect of ~ 12%. 

We account for spatial beam smearing from the PSF 
and velocity broadening from instrumental spectral res- 
olution in our modeling, as we compare this spatially 
and spectrally convolved model disk to the observations. 
Hence, the best-fit quantities represent the intrinsic prop- 
erties derived by properly taking into account these ob- 
servational effects. 

3.2. Fitting technique: an evolutionary approach 

Given the large number of parameters and the poten- 
tial for many local minima in the likelihood/x^ space, 
it is important to use a reliable fitting algorithm, i.e. 
robust and efficient one to find the global extremum. 
We selected for this task genetic algorithms, which are 
heuristic search techniques that provide a means to ef- 
ficiently sample the parameter space, reproducing in a 
computational setting the biologi cal eff ect of evolution 
by natural selection (e.g. Holland. Il975f ). These kinds of 
algorithms have the advantage of being robust and effi- 
cient in finding the global maximum of a given function 
even in a very complex parameter space, with a reason- 
able number of iterations and without the need of a good 
set of initial guesses. It thus represents an ideal tool for 
our specific problem. 

We use a modified version of PIKAIA, a genetic 
algori thm-based optimization routine by Charbonneau 
()1995f l. The algorithm starts with a random population 
of "individuals" that represent a single location in pa- 
rameter space with a complete set of the parameters we 
want to fit. The fitness of each individual is evaluated 
through their x^j the probability to survive in the 
next iteration is set to be proportional to this value. In 
this way the more suitable parameter sets are selected, 
but at the same time the parameter space is efficiently 
sampled by "breeding" between different individuals and 
allowing random "mutations" and "crossovers" between 
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the parameters. Therefore, in each iteration the whole 
population becomes more and more suitable, while still 
sampling the parameter space. When the requested tol- 
erance is reached, or at the end of a given number of 
iterations, the best fitting "individual" is chosen as the 
solution. 

The genetic algorithm allows us to efficiently mini- 
mize the differences between the exponential disk model 
and the observed velocity and velocity dispersion fields 
of the Ha emission line. Simulations on test galaxies 
have shown that leaving the rotation center as a free pa- 
rameters in the fitting significantly worsens the the final 
results, even after increasing the number of maximum 
iterations allowed (see also Shapiro et al. |2008). We 
therefore locate the center using the continuum light dis- 
tribution of the galaxies, which corresponds roughly to 
rest frame R band as observed in the near-IR at 2 ~ 2. 
In the SINFONI data of our galaxies, the continuum is 
detected, although typically at a lower S/N than the 
emission lines. The resulting intensity map is generally 
smoother and distinc t than the emission line (Forster 
Schre iber et al. l2006l Genzel et al. 2008, Shapiro et al. 
120081 ). as it is also the case for local disk galaxies (e.g. 
Daigle et al. |2006). These continuum intensity maps 
show a clear peak in the central region, especially in the 
BzK selected galaxies, and represent a better tracer of 
the mass distribution (dominated by stars) in the galax- 
ies than the more clumpy line emission (tracing active 
sites of star formation). The center is identified by com- 
puting the intensity weighted mean of the light profile of 
the central pixels along both minor and major axes, and 
it is kept fixed in the disk fitting procedure. 

The other parameter that is provided as an input to 
the fitting routine is the scale length of the disk. 
This parameter is in fact better constrained by the 
emission intensity than by the dynamical properties 
of the gas. In this case we used the higher S/N line 
emission image of the galaxies to quantify it. The disk 
Ha distribution has been found to follow the continuum 
for local spiral galaxies, with consistent derived scale 
lengths (e.g. Tacconi et al. J986, Koopmann et al. 2001, 
Hanish et al. 12006} ). From an analysis of different SINS 
and simulated galaxies, we found that the most robust 
way to measure the scale length of the disks given the 
S /N and resolution of the data is to fit a linear Gaussian 
profile along the major axis of the line intensity maps. 
This method has the advantage of using just the bulk of 
the brighter line emission, avoiding the problems related 
to the surface brightness at th e edges of the galaxies. 
Following Bouche et al. ()2007D . we assumed that the 
derived HWHM, once corrected for the observed seeing, 
corresponds to the scale length Rd of the disk. This as- 
sumption is supported by model disk fits, incorporating 
the beam smearing due to the ~ 0.5" FWHM resolution 
of our data. These simulations show that Rd measured 
in this way is likely to be overestimated by no more 
than < 15%. Using available HST/NICMOS continuum 
imaging (see Forster Schreiber et al. in preparation) 
and deep ISAAC imaging in GOODS (data release vl.5: 
|http:/ /www. eso.org/science/goods/releases/20050930/) 
of a subsample of these sources (1 GMASS, 2 K20 
and 5 BX), we find reasonable agreement between our 
measured scale length and the ones derived by fitting 
an exponential profile to these K band images with 



GALFIT (Peng et al. l2002f) . The results are shown in 
Fig.m 

According to simulations, a choice of a scale length 
Rd 15% smaller would provide a dynamical mass 10% 
lower, and a ~ 30% smaller one a mass ^ 20% lower. We 
note that the value obtained for the maximum rotational 
velocity Vmax is instead much less sensitive to a variation 
of the scale length, as the fitting routine has to reproduce 
the observed velocity pattern. However, we also note 
that if a flatter mass distribution is a more appropriate 
description of the structure of these high-z dis ks, e.g . the 
ring- like structures discussed in Genzel et al. (|2008f ). the 
true enclosed masses would then be up to ~ 10 — 30% 
lower than the one obtained with the method presented 
here. 

The inclination of the disk with respect to the line of 
sight I, the position angle on the plane of the sky 9, the 
total enclosed mass Mdyn, the velocity zero point vq at 
the rotation center, and the constant global isotropic dis- 
persion term (T02 (see eq. [2]) , are left as free parameters 
in the minimization. The model disk obtained is then 
rescaled as it would appear at the redshift of the source, 
re-binned to the pixel size of the current observation, 
and smeared using the FWHM of the PSF observed (see 
Sect. [2]). In every iteration of the genetic algorithm, the 
is evaluated for each "individual" from the squared 
difference between the observed velocity and velocity dis- 
persion maps and the modeled ones. The difference in 
each pixel is weighted using the uncertainties on the ve- 
locity and on a derived as described in Sec. [2] at that 
location, e.g.: 



1 



(3) 

where Npi^ is the number of pixels with S/N over the 
threshold (see Sect. 12. 3p used in the analysis. The best 
fitting parameter set is therefore taken as the solution. 
The uncertainties are evaluated following Avni ()1976l ) 
method to calculate confidence limits for a multiple num- 
ber of parameters, studying the variation of around 
the optimized minimum. We assumed as 90% confidence 
interval where the reduced is increased by Ax^ =6.25 
for the three "interesting parameters" (inclination, dy- 
namical mass and (T02) used in the analysis. 

4. RESULTS 

For each galaxy in the sample we show, in Fig. [21 the 
results of the best fit with the exponential disk model, 
for both the velocity and velocity dispersion fields. The 
disk model is able to reproduce the global dynamical pat- 
tern of the galaxies, although in the observed maps small 
deviations from a perfect "spider-diagram" are observed. 
These deviations are clear as residuals in the difference 
image between the observed maps and the model, and in 
the x^ maps, defined as above (eq.[3]). The total reduced 
X^ obtained for the galaxies vary between 0.2 and 20, 
but comparing directly the different x^ values between 
the galaxies is not a reliable criterion for evaluating the 
quality of the modeling for the different sources nor for 
distinguishing ordered motions from disturbed kinemat- 
ics due to interactions. In fact the x^ test is, as expected, 
sensitive to the average S/N of the observation, since the 
value is weighted using the uncertainties on the velocity 
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and dispersion. Since the S/N of the data varies from 
source to source, the separation between rotating disks 
and disturbed kinematics is performed using t he kineme- 
try criteria developed by Shapiro et al. ()2008l . see section 

2)- 

The results of the modeling are summarized in Tab. [21 
The dynamical mass reported is the total mass in the best 
fitting exponential disk at a radius of 10 kpc (~ 1.2" at 
z ~ 2). This radius was chosen to roughly match the 
photometric radius used to measure the "total" fluxes of 
the galaxies in the different filters, which were used in 
the SED fitting to calculate the stellar masses (see fol- 
lowing section) , allowing a meaningful comparison of the 
two quantities. Vmax is the maximum velocity in the ex- 
ponential disk model. The main source of uncertainty is 
the determination of the inclination with respect to the 
plane of the sky. Moreover, the stated errors on the dy- 
namical mass do not include possible systematic effects, 
e.g. due to the value used for Rd or due to the choice of 
the exponential mass profile (see Sect. I3.2| ). 

The measurement of cto2 is mainly constrained from 
the dispersion in the outer part of the disks, where the 
effects due to the beam smearing of the steep central ro- 
tation curve are less prominent. Some of the galaxies 
were too compact, or with insufficient data quality to 
derive a robust measure of this quantity. In one case, 
Q2343-BX389, the irregular dispersion profile was also 
preventing a unique and robust estimate for the isotropic 
dispersion using this technique. In agreeme nt wit h the 
earlier findings of Forster Schreiber et al. (j2006l ) and 
the modeling of Genzel et al. (2003)j the majority of 
the star forming galaxies modeled have a large compo- 
nent of local random motion, with cro2 ^ 30 — 80 and 
a average Vmax/o'o2 — 4.4, with a typical uncertainty 
of ±20 km/s. Low values of the V/sigma ratio in these 
galaxies suggest that a significant fraction of star-forming 
disks at z ~ 2 may be highly turbulent and geometrically 
thick. Concurring evidence from kinematics has now also 
been found in other z ~ 1.5 — 2.5 disk-like sys tems (e.g., 
Wright et al. l2007l van Starkenburg et al. l2008f l and 
is supported by the large thickness 1 kpc) inferred 
from morphological analysis of candidate edge-on disks 
based on high resolution HST imaging (Elmegreen & 
Elmegreen l2006f ). 

The dynamical evidence of the presence of such large, 
rotating and turbulent disks is in agreement with the sce- 
nario in which rapid, smooth gas accretion from the halo 
may play a significant role i n the the formation of mas- 
sive galaxies (Genzel et al. I2008D . This picture is also 
supported by recent broad-band observations of high-z 
populations, which have found a fairly tight relation be- 
tween stellar mass and star formation rate, a sign that 
steady star formation rates, rather than merger induced 
short luminous bursts, domin ate th e stellar assem bly of 
these galaxi es (Noe ske et al. 120071 : Elbaz et al. 120071 : 
Daddi et al. f200 7bD . The most recent dark matte r sim- 
ulations such as the Millennium (Springel et al. !2003) 
show that subsequent mergers are not frequent enough 
to convert all z ^ 2 disks into elliptical galaxies at z = 0, 
and that most of the mass growth is achieved via merg- 
ers less in tense than 10:1 or smooth accretion (Genel et 
al. l2008f l. Recent hydrodynamic simulations (see e.g. 
So mmer- Larsen et al. I2003< . Ocvirk et al. 2008, Dekel et 
al. l2008l ) also show that steady, cold gas streams, which 



penetrate effectively through the shock-heated media of 
dark-matter haloes, are a main accretion mechanism for 
galaxies in this redshift range. 

4.1. Comparison between dynamical and stellar mass 

Thanks to the multiband coverage of the different fields 
from which we drew our SINS targets (see Sect. 12. ip . 
we were able to perform a Spectral Energy Distribution 
(SED) analysis using the available broad-band photome- 
try from the observed optical to near-IR/mid-IR. In par- 
ticular, as described in detail in Forster Schreiber et al. 
(I2008D. we used data from publ ished w orks by Erb et 
al. (12006b', BX), Daddi et al. ('2004a', K20), Kong et 
al. ^006, D3a), Kurk et al. (2006, GMASS), Capak 
et al. ()2007i . COSMOS), as weU as photometry kindly 
provided to us by McCracken et al. (in preparation, 
COSMOS). SED modefing has been carried out by these 
authors but using different model ingredients or assump- 
tions. In order to ensure more uniformly derived prop- 
erties among our sample, we re-modeled all the galaxies 
followi ng th e prescriptions described in Forster Schreiber 
et al. (12004 see also Forster Schreiber et al. 2008). Our 
results do not differ significantly from the published ones. 
We used the solar metallicity Bruzual & Chariot (200;|) 
stellar population synthesis models, with a Chabrier Ini- 
tial Mass Function (IMF ) between 0.1 and 100 Mq, and 
the Calzetti et al. ()1994[ ) reddening law. We considered 
three combinations of star formation history and dust 
content: constant star formation rate (CSF) and dust, an 
exponentially declining star formation rate with e-folding 
timescale of t = 300 Myr and dust (t300), and a dust- 
free single stellar population formed instantaneously. We 
adopted the best of those three cases based on the re- 
duced chi-squared value of the fits (for all the galaxies 
in the sample, this is either CSF or t300). The red- 
shift is kept fixed to the spectroscopic redshift measured 
from the emission lines detected, while the stellar mass 
M,, SFR, star formation history, stellar population age 
and visual extinction Ay are free parameters. We eval- 
uate the uncertainties through Monte Carlo simulations, 
where the input photometry is varied randomly assuming 
the photometric uncertainties are Gaussian. The fitting 
results for the galaxies of the sample are reported in Ta- 
ble ffl 

The reported stellar masses are affected by various un- 
certainties and systematic errors. In particular, different 
stellar population synthesis models do not reproduce a 
consistent picture of evo lution in the rest frame NIR (see 
Marchesini et al. I2008D . For example, Mar aston et al. 
(|2006D suggests that the Bruzual & Chariot l|2003f) stel- 
lar population synthesis models used here may overesti- 
mate the stellar masses up to a factor of '--^ 2 when the 
bulk of the population in the age range between 0.5 
and ~ 2 Gyr. This is due to a different approach in 
the treatment of the thermally-pulsing asymptotic giant 
branch (TP-AGB) phase. However, findings from much 
larger samples, where the effects of SED modeling can 
be better and more reliably investigated on a statistical 
basis, suggest that the correction between the two dif- 
ferent approache s can be negligible or u p to ~ 30% (see 
e.g Ber ta et al. I2007L Marchesini et al. |200S, Wuyts et 
al. l2007l ). Accordingly, using the different library of syn- 
thetic spectra of Chariot & Bruzual (2007), the stellar 
masses estimated for our small sample change individu- 
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ally by no more than 20%, while the median of the whole 
sample remains consistent within 2%. In any case, this 
correction would make the stellar masses lower, and the 
offset in the zero point of the high-z TuUy-Fisher relation 
discussed in Sect. [5] will be correspondingly larger. 

In four objects of the sample there are independent 
indications of the presence of an AGN: dynamical and 
UV spectral evidence suggests th e pre sence of an AGN 
in D3a-15504 (see Genzel et al. 120061 Daddi et al. in 
preparation), BX663 and D3a-7144 show AGN features 
in optical spectroscopy (Erb et al. i2006c and Daddi et al. 
in preparation respectively), while for D3a-6004 a 24^m 
excess was observed by Daddi et al. (in preparation), 
and a large [Nil] /Ha ^ 0.63 ratio was measured in our 
data in the nuclear region of the galaxy (Buschkamp et 
al. in preparation). The inspection of the optical to 
K band available SED for these galaxies suggests that 
this AGN contribution may be biasing the stellar mass 
estimate towards higher stellar m ass in the latter two 
cases (see Forster Schreiber et al. 2008). We therefore 
conservatively consider the stellar mass estimate for D3a- 
6004 and D3a-7144 as an upper limit. 

The stellar masses derived for the sample are between 
l-30x 10^° Mq, with an average of 6.3 10^° Mq confirm- 
ing that all the selected galaxies are drawn from the mas- 
sive star forming galaxy population at high-z. Thanks to 
the SED results and our dynamical model, we are able 
to compare the stellar and dynamical mass of the galax- 
ies. In Fig. [3] these are compared for galaxies in which 
multi-band photometry is available to constrain the SED 
fitting. We observe remarkable correlation between the 
two quantities: the Spearman correlation test gives a 
probability P — 2 ■ 10^^ that the two masses are uncor- 
related. 

The resulting average dynamical mass at i? = lOkpc 
is ~ 3.3 times larger than the average stellar mass, ex- 
cluding the sources with a possible AGN contribution. 
Therefore, given a dark matter contribution of ~ 40% at 
this radius t o the t otal dynamical mass of such disks (see 
Genzel et al. l2008l ). we infer a large fraction of gas in the 
sample galaxies {Mgas M^, or ^ 30%), although with 
very large uncertainties. This estimate may be compared 
with a gas mass derived from the star formation rate sur- 
face density u sing t he Schmidt-Kennicutt relation (see 
Bouche et al. I2007D . In this way we derive an average 
gas contribution of 23% to the total dynamical mass of 
the galaxies, where the difference can be accounted for 
by the large uncertainties in both our estimate and the 
extinction correction needed to derive the star formation 
rates from the detected Ha fluxes. The gas fraction esti- 
mat e is al so in agreement with the findings of Daddi et 
al. ((2008'), who directly observed molecular gas in two 
z ^ 1.5 BzK galaxies, deriving gas fraction as high as 
Mgas ~ M^,. An alternative possibility is that the stel- 
lar mass is significantly underestimated for most of the 
sample, presumably because of a faded old stellar pop- 
ulation that is not detected in the multicolor SED, but 
would still contribute to the total dynamical and stellar 
mass (see the di scussi on by, e.g. Shapley et al. 1200 iL 
Papovich et al. l2001l and for the galaxies with avail- 
able IRAC photometry Wuyts et al. 20Q3), or that the 
dynamical mass is overestimated due to the systematic 
effects discussed in Sect. 13.21 

Although the stellar population age is not strongly con- 



strained by the SED fitting due to the uncertainty on the 
star formation history and the degeneracy with extinc- 
tion, we note that the objects with the higher Mdyn/M,, 
generally correspond to younger best fitti ng ages , sim- 
ilar to what was discussed by Erb et al. (2006b), and 
that none of the youngest galaxies has a small mass ratio 
(see Fig. 0]). The significance of the correlation accord- 
ing to the Spearman test gives a probability P = 0.08 
that the two quantities are uncorrelated. The scatter is 
in fact high, probably suggesting that these galaxies are 
not closed systems, and that an important contribution 
from continuous gas accretion must be considered. This 
again suggests that the galaxies in the sample are star 
forming disks with significant gas fraction and gas ac- 
cretion, where the gas consumption correlates with the 
stellar population age. 

5. THE Z ~ 2.2 TULLY-FISHER RELATION 

The relation between luminosity (or stellar mass) and 
maximum rotational velocity of di sk galaxies has long 
been known (TuUy & Fisher [1977|) . This TuUy-Fisher 
relation (TFR) has been used as a dista nce indic ator fo r 
disk galaxies (e.g. TuUy & Pierce l2000l, SandagelUll), 
but also as a key feature for understanding the structure 
and evolution of these galaxies, as it links directly the 
angular momentum of the dark halo with the luminosity 
(or mass) of the stellar population in the disk. In fact, ac- 
cording to the stan dard model for disk galaxy formation 
(Fall & Efstathiou Il980j ), disks form out of gas cooling 
down from a hot halo associated with the dark matter 
potential well, maintaining its specific angular momen- 
tum a nd se ttling in a rotationally supported disk (Mo 
et al. |1998[ ). Consequently, the structure and dynam- 
ics of disk galaxies are expected to be closely correlated 
with the properties of the dark matter halo in which they 
are embedded, and the observable structural parameters 
of disks at different cosmic epochs can be used as trac- 
ers of the properties of the corresponding halos. The 
evolution of the stellar mass TFR is expected to be re- 
lated both to the conversion of gas into stars and to the 
inside-out growth of the dark matter halo by accretion. 
In fact, while the extent of the halo grows considerably 
with time, the circular velocity of the halo grows less 
keeping the rotation curve approximately flat to larger 
and larger distances. The accretion of the dark matter 
is followed by accretion of baryonic gas, which is subse- 
quently converted into stars by ordinary star-formation 
in the disk. However, the details of the process and the 
amount of evolution expected depends strongly on the 
model assumptions, accretion mechanism adopted and 
the timescale needed to convert the gas into stars. Any 
successful model of disk formation should then be able 
to reproduce the slope, zero-point, scatter and redshift 
evolution of the TFR and other scaling relations. How- 
ever, reproducing simultaneously the TFR, the shape 
and normalization of the luminosity function, and the 
observed sizes, metallicities and colors of disk galaxies 
is still a ch allenging task for current models (Gourteau 
et al. l2007l ). Therefore, as the normalization, slope and 
evolution of the TFR depend strongly on the prescrip- 
tion used for the star formation and feedback, as well as 
on the disk formation mechanism and cosmological pa- 
rameters, it represents a cruc ial test for galaxy evo lution 
models (e.g. Silk et al. Il997t Steinmetz & Navarro ll999l : 
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van den Bosc h l2000l l2002t Sommer-Larsen et al. 120031 : 
Button et al. l2007f ). 

In recent years an increasing number of dynamical ob- 
servations of disk galaxies in the local and intermediate 
redshift Universe have allowed us to begin exploring the 
evolution of the TFR with redshift. There are still dis- 
crepant results on a possible evolution at intermediate 
redshifts of the tight relation obser ved in local samples 
of dis k galaxies (e.g. Haynes et al. 119991 Pizagno et al. 
12007^ between the absolute magnitude and the maximum 
rotational velocity. For example, while Vogt et al. (,1996) 
reported very little evolution of the B-band TFR up to 
z ~ 1, other groups (see e.g. Simard & Pritchet I1998L 
Ba rden e t al. 2003, Nakamura et al. |2006| and Bohm et 
al. 12004') found a strong brightening of ~ 1 — 2 mag in 
B-band luminosity over the same redshift range. 

Furthermore, the interpretation of a possible evolution 
of a luminosity based TFR is difficult as both the lumi- 
nosity and the angular momentum might be evolving at 
the same time. For this reason, the stellar mass TFR, 
which correlates the stellar mass and the maximum rota- 
tional velocities of disks, offers a more physically robust 
comparison as i t invo lves more fundamental quantities 
(B ell fc de Jong|2001|, McGaugh et al. l2005l Pizagno et 
al. 12005!). In this context, Flores et al. (2006'), Kassin et 
al (2007j) and Conselice et al. ([2005) have found no evo- 
lution in both the slope and zero point of the stella r mass 
TFR up to z ~ 1. On the other hand, Puech et al. (|2008D 
detect an evolution of the stellar mass TFR zero point of 
0.36 dex between z ^ 0.6 and z = 0, using a sample of 18 
disk-like galaxies observed with the integral field spectro- 
graph GIRAFFE at the VLT. The di fferences with the re- 
sults obtained by Flores et al. (j2006l ) on a similar dataset 
are attributed to both a difference in the reference local 
relation and a more accurate measurement of the rota- 
tion velocity at high redshift. A very mild evolution was 
also found in cosmological N-body/hydrodynamical sim- 
ulatio ns of d isc galaxy evolution (Portinari & Sommer- 
Larsen [2007l see Sect. 15. l( ). where the galaxy growth with 
time since z ~ 1 is happening mostly along the TFR. 
The growth of the disk baryonic mass by infall of gas 
and mergers between z ~ — 1 is therefore observed to 
be accompanied by a correspondingly large increase of 
the rotation velocity, so that the evolution of individual 
objects occurs mainly along the TFR. 

5.1. The observed z ~ 2 Tully-Fisher relation 

With the available modeled kinematics of z ~ 2 galax- 
ies from our SINS data, we are now able to push, for 
the first time with a sizeable sample, the investigation of 
the evolution of the TFR to even higher redshift. When 
comparing the evolution of the high-z TFR with that 
inferred from local (or lower redshift) samples of star 
forming galaxies, one must be aware that we are compar- 
ing different classes of objects, which are not necessarily 
linked from an evolutionary point of view. Therefore any 
evolution in the relation can not be interpreted straight- 
forwardly as evolution of any individual galaxy, as from 
z — 2 to z — part of our galaxies may undergo signifi- 
cant (morphological/dynamical) evolution driven by effi- 
cient secular/internal dynamical processes and/or merg- 
ing (cf. Genzel et al. 2008, Gencl ct al. 2008). Rather, 
evolution of TFR represents the evolution of ensemble 
properties of disk-like populations at different redshifts. 



The available full 3D coverage of the dynamical fea- 
tures will not only provide a more complete and detailed 
view of the kinematics, but will also allow us to avoid 
the biases and limitations due to slit spectroscopy, such 
as misalignment of the slit and the major axis, and to 
identify disturbed rotational patterns avoiding the con- 
tamination of our disk-like sample with more complex 
dynamics. Using similar 3D SINFONI data, Starken- 
burg et al. (|2008f l have presented the position of a single 
galaxy at z = 2.03 on the stellar mass TFR, yielding an 
offset of 2a) from the local one (see Fig. [5]) . In view 
of the uncertainties of the quantities involved, which are 
not insignificant at high-z, larger samples are required to 
constrain more robustly any offset in the TFR at z ~ 2. 

We have therefore used the best fitting Vmax from our 
dynamical modeling, together with the stellar mass from 
SED fitting discussed in sect. 14.11 to build the TFR for 
our sample, shown in Fig. [5l The solid line sho ws th e 
z = best fit relation obtained by Bell & de Jong (|2001h . 
Both are corrected to convert their "diet" Salpeter and 
Salpeter IMF to the Chabrier IMF used here, using a fac- 
tor 1.1 9 and 1.7 respectively. We note that Pizagno et al. 
()2005l ) have obtained a shallower (M^ oc T/3 .05-) jQ^g^j g^gj_ 
lar mass TFR than Bell & de Jong ()2001l . (xV^-^), 
but the two relations agree at V — 200 km s~^. How- 
ever, they used the circular velocity at 2.2 Rd instead of 
Vmax to estimate of the rotational velocity, and differ- 
ences in the sample selection ma y also account for the 
different result. McGaugh et al. (|2005| ) obtained an in- 
termediate slope, Mbar oc but they explored the 
total baryonic TFR, including the contribution of the 
gaseous mass component. In this case, higher gas frac- 
tion in lower mass gala xies m ay explain the flatter slope. 
Recently Meyer et al. (l2008l ) presented the local stellar 
mass TFR for a large Hl-selected sample of galaxies, ob- 
taining a slope consistent to the Bell & de Jong (2001]) 
value adopted here for comparison with previous studies. 

The circles in Fig. [5] show the data for the z ~ 2.2 
galaxies from SINS, while the filled triangles the z ^ 1.5 
sources. The open circle is the z 2 galaxy presented 
m van Starkenburg et al. (|2008D . The average uncer- 
tainty is shown as an error bar in the lower right cor- 
ner. This does not include possible systematic effects 
both in the SED fitting (see Sect. 14. ip and in the dy- 
namical modeling (see Sect. The z ~ 2.2 galaxies 
show an evident evolution in the zero point of the TFR. 
The dashed line shows the best fitting relation keeping 
the sl ope fi xed at the z = value found by Bell & de 
Jong (|2001h . The current data are in fact not yet able 
to fully constrain the slope of the relation due to the 
limited statistics and mass range. The best fitting zero 
point is ZP2.2 = —0.09 ± 0.11, where the error is the Icr 
uncertainty, to be compared with the local ZPq — 0.32, 
giving a significance of ~ 3.7a. The z = 2.2 relation can 
be written as log{M^) = —0.09-1-4.5 x \og{Vmax)-, finding 
an offset in log{M^,) of 0.41 for a given rotational veloc- 
ity. We do not find evidence for significant evolution in 
our z ~ 1.5 sample with respect to the local relation, al- 
though the data are still too sparse in this redshift range 
to place a firm constraint. The evolution observed at 
z ~ 2.2 is larger than the intermediate 0.36 dex found 
at z ~ 0.6 by Puech et al. (2008), although they use a 
different local reference, suggesting an increasing of the 
evolution of the TFR with redshift. 
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5.2. Comparison with galaxy evolution models 

Interestingly, the amount of observed evolution in our 
sample is consistent with the factor ^ 1.25 of increased 
velocit y for a given stellar mass predicted by Somerville 
et al. ( 20081 ) in their semi-analytic formation model of a 



disk with M, = 6 x 10^° Mq at this redshift. They follow 
the classical paradigm for disk formation within massive 
extended d ark m atter halos initially propose d by Fall & 
Efstathiou (|1980( ) but use a Navarro et al. ()1997l ) mass 
density profile instead of a singular isothermal sphere as 
halo mass density profile, account for the expected evo- 
lution wit h tim e of the profile concentration Cmr (Bul- 
lock et al. l2001f ). and for the effect of self-gravity of the 
baryons in the central part of the halo (so called "adi- 
abatic contraction"). They find that with this refined 
model the disk rotation velocity Vmax at fixed disk mass 
evolves much more gradually than the halo virial velocity 
at fixed halo mass, remaining nearly constant and show- 
ing the largest evolution between z ~ 1.5 — 3. If this 
prediction of the model is correct, our sample therefore 
probes the best suited cosmic epoch to test the evolution 
of the TFR. 

We also compare the observed TuUy-Fisher relation 
with the Smoothed Particle Hydrodynamic (SPH) sim- 
ulations of disk formation and evolution in a cold 
dark matter hierarchical clustering scenario by Sommer- 
Larsen et al. (in preparation) and Portinari & Sommer- 
Larsen (j2007t l. They include in their modeling star for- 
mation and stellar feedback effects to delay the collapse 
of protogalactic gas clouds and allow the infalling gas to 
preserve a large fraction of its angular momentum. Pre- 
vious numerical simulations have in fact shown that when 
only cooling processes are included, the infalling gas loses 
too much angular momentum, resulting in disks that are 
much smaller than the obser ved o nes (Navarro & White 
1 1994 Navarro & Steinmetz l2000( l. Moreover, chemical 
evolution with non- instantaneous recycling, metallicity- 
dependent radiative cooling and the effects of a meta- 
galactic ultraviolet field, including simplified radiative 
transfer, were implemented in the simulations. This is 
an improved version of th e models developed by Sommer- 
Larsen et al. ()1999Ll2003h to address the angular momen- 
tum problem in disk formation models and provide realis- 
tic z = galaxies. It has proven successful in reproducing 
better than previous numerical models the observed sizes 
of disk galaxies, the luminosity profile of both disks and 
bulges, bulge to disk ratios, integrated colors and TFR 
of local galaxies, as well as the observed peak in the cos- 
mic star formation rate at z ~ 2 (Sommer -Larsen et 
al. l2003f l. The models are characterized by a formation 
mechanism where a large fraction of the gas is accreted 
onto the forming disk through both rapid cold gas ac- 
cretion in filamentary structures (dominant at z > 2) 
and gradual cooling flows from a surrounding hot phase 
(dominant for z < 2), without abrupt changes of angular 
momentum orientation, rather than by rapid mergers of 
massive cold clumps. These simulated galaxies at z ~ 2.2 
are shown as open diamonds in Fig. O Again, we obtain 
very good agreement between our data and the predicted 
position of z ~ 2.2 galaxies on the TFR for a range of 
stellar masses. In the model the predicted zero point off- 
set is produced by a relevant inside-out mass accretion 
between z ~ 2.2 and z ^ 1, that do not correspond to 



an equal increase of the rotation velocity of the galaxies. 
This issue will be investigated in a forthcoming paper 
(Sommer-Larsen et al. in preparation). However, if the 
TFR evolution is driven mainly by dark matter as sug- 
gested by the models, this may suggest that the galaxies 
we are observing at high-z are "stable" , in contradiction 
with inferences from our detailed studies of several of 
our disks which appear to be glo bally unstable to star 
formation (see Genzel et al. l2008|). 

5.3. The S'o.5 Tully-Fisher relation 
Kassin et al. (|2007t l have constructed a stellar mass 
TFR at z < 1 using the S'o.5 = ^0.5 • V^^^ + erg esti- 
mator, in order to take into account disordered or non- 
circular motions through the gas. In this way, they found 
a tighter relation up to z = 1.2, with no significant evolu- 
tion, and conclude that the observed scatter in the TFR 
is partly due to disturbed and peculiar morphologies, cor- 
responding to lower rotation velocities for their masses. 
In Fig. [S] the result obtained for our z '--^ 2 sample us- 
ing the 5'o.5 index is shown, using the fitted CT02 as a 
measure of the intrinsic dispersion of the gas not due 
to rotation. Although the a used here does not exactly 
match the definition adopted by Kassin et al. (2007), 
it is the best indicator of turbulent and disturbed mo- 
tion obtained from our modeling. Using the S'o.5 index 
we also detect an evolution of the relation between z ~ 1 
and z 2.2, further supporting the claim of an evolution 
of the TFR. However, due to the different definition of a 
used here, we do not attempt to quantify this evolution. 
The resulting scatter is slightly smaller than the already 
remarkably low one obtained for the classical TFR (see 
Fig. [5]). We suggest that this is due to the optimal se- 
lection of our rotational supported systems at z ~ 2, as 
the SINFONI 3D data are ideal to ensure a much lower 
contamination from complex dynamics and irregular mo- 
tions than in longslit datasets. 

6. CONCLUSIONS 

We have presented the dynamical modeling of the full 
2D velocity and velocity dispersion fields of the Ha emis- 
sion line for 18 star forming galaxies at z ~ 2, ob- 
served with the integral field spectrometer SINFONI in 
the near-IR. The galaxies are selected from the larger 
sample of the SINS survey, which probes galaxy popula- 
tions responsible for the bulk of the cosmic star forma- 
tion activity and stellar mass density at these redshifts, 
based on the prominence of ordered rotational motions 
with respect to more complex merger induced dynamics. 
The quality of the data allows us to model the observed 
dynamics with suitable simulated exponential disks, and 
derive the total dynamical mass and maximum rotational 
velocity of these galaxies. 

We find a good correlation between the dynamical 
mass and the stellar mass derived through SED fitting, 
deriving that large gas fractions {Mgas ~ A/*) are re- 
quired to explain the difference between the two quanti- 
ties. We also find that the galaxies with younger stellar 
population age have genera lly larger gas fractions, as re- 
ported by Erb et al. (|2006br ). although the large observed 
spread could be related to a significant contribution of 
continuous gas accretion. 

We used the derived stellar mass and the Vmax of the 
best fitting disk model of each galaxy to construct for 
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the first time a TuUy-Fisher relation at z ^ 2.2. This is 
the redshift range where the largest evolution may be ex- 
pected according to some recent simulations (Somerville 
et al. l2008f) . The full 2D spatial coverage of our data 
provides a view of the kinematics in unprecedented de- 
tail at these redshifts, with which we can avoid the bi- 
ases and limitations of slit spectroscopy in studying the 
dynamics of high-z sources. Moreover, the full 2D dy- 
namical information allows us to select a sample with 
a lower contamination from irregular and complex dy- 
namics, which is better suited for studying the evolu- 
tion of the TFR in rotationally dominated systems. We 
find a remarkably low scatter in the TuUy-Fisher rela- 
tion in our sample, and an evolution of the zero point 
in log{M^,) of 0.41 ± 0.11 for a given rotational veloc- 
ity, which is fully consistent with the predictions of the 
latest N-body/hydrodynamical simulations of disks for- 
mation and evolution. To reproduce the observed prop- 
erties of the local and high-z TFR, the SPH model used 
as comparison (Sommer-Larsen et al. il9_99!, l2003l and 
in preparation) invokes a formation mechanism where 
a large fraction of the gas is accreted onto the forming 
disk through both rapid cold gas accretion in filamen- 
tary structures and gradually in a cooling flow from a 
surrounding hot phase, rather than by rapid mergers of 
massive cold clumps. This is in qualitative agreement 
with recent results showing that smooth gas accretion 
from the halo is required to reproduce the observed prop- 
erties of a large fractio n of th e z ~ 2 galax ies (see e.g . 
Forster Schreiber et al. l2006l : Genzel et al. l2006l 120081 



Shap iro et al. 120081: No eske et al. l2007t Elbaz et al. 
I2007t Daddi et al. l2007bl ). As aheady suggested by the 
velocity-size relation obtained for a larger sam ple of SINS 
UV/optically selected galaxies (Bouche et al. l2007l ). the 
angular momentum properties of these z ~ 2 galaxies 
are very similar to local spirals and intermediate redshift 
disks. However, these z ~ 2 disks show a much larger 
component of random local gas motion {V/ao2 ~ 4.4, 
and up to an order of magnitude less than z = spi- 
rals), implying high turbulence and scale heights. The 
random motions could be the result of feedback from the 
intense star formation (~ 50 — 200 Mq /yr) , or be driven 
by the accreting gas as it settles onto the forming disk 
(see Genzel et al. |2008[) . 

The results obtained confirm that the capabilities of 
Integral Field Spectroscopy can provide unique insights 
in the formation and evolution of high-z galaxies, and 
that large, gas rich, rotating disks, represent a significant 
population in the z ~ 2 universe. 
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Fig. 1. — Comparison between the scale length measured from the line emission in our SINFONI data with the value obtained by 
fitting a Sersic profile to HST and ISAAC K-band continuum images available for a subsample of galaxies (5 BX, 2 K20 and 1 GMASS). 
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Fig. 2. — Kinematic fitting of the galaxies of the sample. For both the velocity (upper panel) and velocity dispersion (lower panel), 
we show from left to right the observed maps, the best fitting model maps and their differences. It can be seen how the global rotational 
pattern is typically well reproduced, although local deviations from the ideal model are highlighted in the difference. The Ha line map is 
shown in the right panel for each galaxy, with the intensity scale normalized to the brightest pixel. 
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(continued) - Kinematic fitting of the galaxies of the sample. 
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Fig. 2. — (continued) - Kinematic fitting of the galaxies of the sample. 
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(continued) - Kinematic fitting of the galaxies of the sample. 
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Fig. 2. — (continued) - Kinematic fitting of the galaxies of the sample. 
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9.5 10 10.5 11 11.5 12 

Fig. 3. — Comparison between dynamical and stellar mass for the galaxies of the sample with available multiband photometry. The 
dynamical mass is the total mass in the best fitting model up to i? = 10 kpc, while the stellar mass is derived by SED modeling assuming a 
Chabrier IMF. The filled triangles are the z ~ 1.5 galaxies from our sample, while the circles are at 2 2.2. The solid line indicates equal 
mass, while the dashed one indicates a dynamical mass 3.3 times larger. The error bar in the lower right corner represents the average 
fitting uncertainties in both directions. 
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Fic;. 4. — The mass ratio M^y^^/Mf is plotted against the best fitting age from SED modeling. Tlie filled triangles arc the z ~ 1.5 
galaxies from our sample, while the circles are at z ~ 2.2. Galaxies with younger stellar populations show a higher Miiyj^/Mf, suggesting 
larger gas fractions. However, the large observed scatter may suggest a significant contribution from continuous gas accretion. The error 
bar in the lower right corner represent the average fitting uncertainties in both directions. 



Modeling of the dynamics of z ~ 2 galaxies 
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loe(V ) 

&^ max' 

Fig. 5. — The stellar mass TuUy-Fisher relation ai z ^ 2. The filled triangles are the z ~ 1.5 galaxies from our sample, while the filled 
circles are at 2; ~ 2.2. The empty circle show the z = 2.03 disk galaxy F257 observed by van Starkenburg et al. (2008). The error bars 
in the lower right corner represent the average fitting uncertainties of the model maximum velocity and of the stellar mass from the SED 
fitting, not accounting for systematic errors. The solid line is the 2 = relation from Bell & de Jong (j2001|). corrected to a Chabricr IMF. 
The dashed line is instead the best fitting zero point to the z ^ 2.2 observed galaxies, keeping the slope fixed at the z ~ value from 
Bell & de Jong (j20bl1 '). The open diamonds are the 2 ~ 2.2 simulated galaxies of Sommer-Larsen et al. (2008), where the bars show the 
variation in the rotational velocities in the model assuming a 2 2.2 scale length ranging from the 2 = value to half of it. 
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Fig. 6. — The stellar mass TuUy-Fisher relation at 2 ~ 2 using the 5o.5 = vTO-S ■ ^rot + "^0) estimator introduced by Kassin et al. II2007I) , 
in order to take into account disordered or non-circular motions through the gas. The filled triangle s are t he 2 ~ 1.5 galaxies from our 
sample, while the circles are at z ^ 2.2. The dashed line is the 2 ~ 1 relation derived by Kassin et al. 1 120071) . 
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TABLE 1 

The galaxy sample selected for the dynamical modeling and their properties, z is the redshift derived from the Hq 

EMISSION line, mx THE APPARENT K BAND MAGNITUDE, Select. SHOWS THE CRITERIA USED TO SELECT THE GALAXY (SEE SeCT. [TTI l. Scale 
AND Exp. time THE SINFONI pixel scale and integration time used for the OBSERVATIONS, M* AND Age ARE THE STELLAR MASS 
AND AGE OF THE STELLAR POPULATION DERIVED FROM THE SED FITTING OF THE SOURCES (SEE SeCT. [ITTI) . Ha IS THE TOTAL LINE FLUX 
MEASURED IN THE SINFONI DATA AND SFR IS THE STAR FORMATION RATE DERIVED ACCOR DING TO KeNNICUTT ET AL. (loos') FROM THE 
Ha FLUX, CORRECTED FOR A FIXED EXTINCTION OF Av = 0.8 USING A CALZE TTI ET AL. 1110941) LAW AND FOR A ChABRIER IMF (SEE 

FORSTER SCHREIBER ET AL. 120081 FOR DETAILS). 
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TABLE 2 

The RESULTS OF THE DYNAMICAL FITTING OF THE GALAXIES USING THE GENETIC MINIMIZATION, i IS THE DERIVED INCLINATION OF THE 
GALAXY WITH RESPECT TO THE PLANE OF THE SKY; Mjy„{10 kpc) IS THE DYNAMICAL MASS AT 10 KPC (~ 1.2"), WHILE M^y„{R^) AT THE 
SCALE LENGTH R^] (Tq IS THE ISOTROPIC VELOCITY DISPERSION; ADDED THROUGHOUT THE DISK (SEE TEXT. SOME GALAXIES WERE TOO 
COMPACT, OR WITH DATA QUALITY NOT GOOD ENOUGH, TO DERIVE A ROBUST MEASURE OF THIS QUANTITY); R^ THE DISK SCALE LENGTH; 
Vmax THE MAXIMUM VELOCITY IN THE BEST FITTING DISK MODEL. THE RELATIVE UNCERTAINTIES ON Mdyn(Rd) ARE THE SAME AS ON 
A^dl/n(10 kpc), WHILE THE ONES ON Vmax ARE 1/2 THOSE ON M^y„{W kpc). 



